Synchronous and Asynchronous Recursive Random Scale-Free Nets 
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We investigate the differences between scale-free recursive nets constructed by a synchronous, 
deterministic updating rule (e.g., Apollonian nets), versus an asynchronous, random sequential 
updating rule (e.g., random Apollonian nets). We show that the dramatic discrepancies observed 
recently for the degree exponent in these two cases result from a biased choice of the units to be 
updated sequentially in the asynchronous version. 
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Stochastic scale-free graphs have been found in a host 
of natural and manmade phenomena (the internet and 
the World Wide Web, networks of flight connections, 
of social contact, of predator prey, of metabolic re- 
actions, eta] and have attracted much recent atten- 
tion P, 0, IE EJ IE El ln this context, deterministic 
scale-free graphs have been very useful as exactly solv- 
able models of the stochastic scale-free nets encountered 
in everyday life 0,0, El El EI- Most deterministic scale- 
free nets are recursive: they are constructed by repeated 
iteration of a fixed set of rules. It is, indeed, their recur- 
sive character that makes them particularly amenable to 
analysis. 

Recently, studies of random variations of deterministic 
nets have yielded interesting results. Consider, for exam- 
ple, the deterministic Apollonian net. Its construction 
begins with 3 interconnected nodes (a complete graph 
of order 3, K3) and at each iteration a new node is 
connected to the vertices of every existing 3-clique (a 
K 3 subset), omitting those 3-cliques that had already 
been updated in previous iterations. In the random ver- 
sion the 3-cliques to be updated are selected at random, 
one at a time. In both cases there results a scale- free 
graph, but the deterministic net has degree exponent 
7 = 1 + In3/ln2 0, El 13 E]L as opposed to 7 = 3 
of the random construct |LJ, Il5| . Similar discrepancies 
are observed for the deterministic network of Dorogovt- 
sev, Goltsev, and Mendes Q (where at each iteration a 
node is connected to the endpoints of each existing link) 
versus the random version (nodes are connected to the 
endpoints of randomly selected links, one at a time), as 
well as for other types of recursive graphs ^(| . Such dra- 
matic differences between synchronous (or parallel) ver- 
sus asynchronous (or random sequential) updating are 
surprising. Barring accidental symmetries, in dynamical 
lattice models and interacting particle systems the two 



kinds of updating usually yield the same kind of behav- 
ior and critical exponents |l7| . 

In this paper we argue that the reason behind the dra- 
matic differences between synchronous and asynchronous 
updating stem from a biased choice of the units to be 
updated in the sequential constructions. We focus on 
a particularly simple scale-free tree whose deterministic 
construction involves connecting new nodes to the end- 
points of each of its links. We then consider two meth- 
ods of sequential updating that differ from each other 
only in the selection rule of the links to be updated. 
The degree exponents obtained in these two ways differ 
from each other (as well as from that of the deterministic 
tree), demonstrating our point. A third method of asyn- 
chronous updating is presented, that avoids the pitfall of 
biased selection and yields the same degree exponent as 
for synchronous updating. We briefly explore the differ- 
ences between random and deterministic constructs even 
when the degree exponents obtained in the two methods 
agree with one another. 
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Synchronous Updating 

Consider the deterministic tree of Fig. QJi, obtained by 
the following procedure : Starting from K 2 at genera- 
tion n = 0, construct successive generations by attaching 
nodes of degree one to the endpoints of each existing 
link (Fig. [TJ>). The tree that emerges in generation n 
has two hubs (the nodes of highest degree) of degree 2™. 
An alternative way for constructing the tree consists of 
doubling the degree of each existing node, from k to 2fc, 
by attaching to it k single-degree nodes (Fig. QJ;). Yet a 
third method, which highlights the self-similarity of the 
tree, consists of producing 3 replicas of generation n and 
joining them at the hubs (Fig-QJl). 

It is clear that all nodes have degrees that are powers 
of 2. Let N n (m) be the number of nodes of degree 2 m in 
generation n. Let N n = J2 m -^«.( m ) be the total number 
of nodes (the order) in generation n. Let M n be the 
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FIG. 1: Synchronous recursive scale- free tree and methods 
of construction, (a) Generations n = 0,1,2,3 of the tree. 

(b) First method of construction: to each of the endpoints 
of every link in generation n connect a node of degree one. 

(c) Second method: to each node of degree k in generation n 
add k new nodes of degree one. (d) Third method: to obtain 
generation n+ 1, join three copies of generation n at the hubs 
(the nodes of highest degree) . 



number of links (the size) in generation n. We have, 

M n = 3M„_x , Mo = 1 , (1) 

(seen most easily by the first method of construction), 
from which follows that 



Mn = 3" . 

Since the graph is a tree, 

N n = M n + 1 = 3" + 1 . 

Also, 

N n (m) = N n -x{m - 1) + 2 • T^S^o 
leading to 



(2) 



(3) 



(4) 
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(5) 



This corresponds to a scale-free degree distribution of 
degree exponent 7 = 1 + In 3/ In 2. 



Asynchronous Updating 

Let us now explore the consequences of asynchronous, 
or random sequential updating, and how it differs from 
synchronous updating. Based on the first method of con- 
struction (Fig-IIJ)), at each time step, t, we choose a link, 
randomly, and connect a new node to each of its end- 
points. We intend to show that the differences between 
synchronous and sequential updating arise because of bi- 
ases in the selection rule of the link to be updated. To 
this end we consider the following two rules: (a) Select 
one of the M(t) links in the net, randomly, with equal 
probability, or (b) Select a node (among the N(t) nodes 
of the net), randomly, then select one of its neighbors, at 
random, and pick the link that connects the two nodes. 
A detailed analysis of these two rules and how they differ 
in the selection of specific links is given in Appendix lAl 

Consider how the degree of node i, fci(t), changes with 
time, by method (a). Since each of the ki links lead- 
ing to node i is selected with probability 1/M(t), the 
probability that ki — > fci + 1 in the next time step is 
ki(t)/M(t). At each time step we add two links to the 
tree, so M(t) = 2t — 1 (we begin with a single link at 
time step t = 1). Hence, in the long time asymptotic 
limit changes in ki are given by 



dki 
~dt 



2f 



(6) 



The initial condition for this equation is ki(ti) = 1, that 
is, we assume that the node was introduced to the tree at 
time i, as a node of degree one (like all newly introduced 
nodes). The solution is then 



(7) 



It follows from Q that the probability that ki is larger 
than k, is 



(8) 



However, since node i could be introduced in any of the t 
steps with equal probability, the probability that tj, < T 
is T/t, so that = 1/fc 2 . The degree distribution for 
large k then follows: 



Pa(k) = -^X(k) - fc- 3 



(9) 



We conclude that the random sequential construction 
of the tree by method (a) leads to a scale-free degree 
distribution of degree exponent 7 = 3, different from 
7 = 1 + In 3/ In 2 of the deterministic tree. 

But what if we select the links by method (b)? In 
this case node i may be the first of the two nodes to 
be selected in step t. This may happen with probability 
1/N(t). Node i may also be the second node selected. 
The probability that we reach i through a randomly se- 
lected node is ki/N(k). The degree of node i increases 
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from hi to ki + 1 regardless of whether it is picked first or 
second. That is, the rate of increase is (l/N)(l + ki/ (k)). 
But N(t) = 2t, while (jfe) = 2M/N = {At - 2)/(2i) -> 2, 
as t — > oo. It follows that in the long time asymptotic 
limit 



dki 



1 

2t 



(10) 



From here we proceed exactly as for method (a), this 
time obtaining 



P b (k) 



(11) 



for large k. Once again the degree distribution is scale- 
free. The degree exponent 7 = 5 is not only different 
from that of the deterministic tree but also differs from 
7 = 3 obtained by method (a). 
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FIG. 2: Random sequential trees constructed by selecting 
links at random (x), or by picking the link between a ran- 
domly selected node and a neighbor (o). The trees consist of 
10 7 nodes. The straight lines of slope 3 and 5 are shown for 
comparison. 

In Fig. |21 we show the degree distribution of random 
sequential trees constructed by the two methods. The 
simulation data are consistent with the degree exponents 
predicted above, though one could argue that larger sim- 
ulations are needed for the distributions to converge to 
their long time asymptotic limit 19] . At any rate, the 
simulations demonstrate our point that a bias in the se- 
lection of the updated units is responsible for the differ- 
ences between different kinds of random sequential nets, 
and between random sequential nets and synchronous 
nets. The differences from the deterministic tree, in our 
case, result from the fact that picking links at random (by 
either of the two methods considered here) favors links 
with endpoints of higher degree. 

A third method for building the random sequential tree 
is based on the recursive technique of Fig. where at 
each iteration the degree of all existing nodes is doubled 
by attaching to them new nodes of degree one. In the 
random sequential version the degree of only one ran- 
domly selected node is doubled at each time step. The 



node to be updated is picked with equal probability from 
among the N(t) nodes of the net. We now show that 
a random net constructed by this method achieves the 
same degree exponent as the deterministic net. 

We cannot resort to the same technique employed for 
the previous two random sequential constructs because 
now ki — ► 2ki in a single step, and the increase is too 
large to allow for a continuous approximation (we are 
trying to study fc, large). Instead, we resort to discrete 
rate equations. 

The degrees in the random tree are powers of 2, just 
as in the deterministic tree. Let N t (m) be the number 
of nodes of degree 2 m at time t, and let N t = J2 m Nt(m) 
be the order of the tree, then 



N t+1 (m) = N t (m) + 



N t +i(P) =N t (0) 



N t (m-l) N t (m) 
N~ t ~N t : 

,N t (m) 



m > 1, 



N t (0) 
N t 



(12) 



Nt 



The long time asymptotic limit can be obtained by 
making the ansatz: N t (m) — > b m t and Nt — * Bt, as 
t — > 00 (B = X) m k m ). Substituting in the first of 
Eqs. I|12fl . we learn that b m = & m _i/(l + B), leading 
to b m = &o/(l + B) m . The second equation tells us that 
B = 2, and the remaining 60 is obtained from the con- 
straint J2 m b "i = B. We get 



N t (m) 



}m+l 



t 



(13) 



This is essentially the same distribution as (JjJJl, corre- 
sponding to 7 = 1 + In 3/ In 2. This shows that there 
exist ways of choosing the units to be updated such that 
the differences between random sequential nets and de- 
terministic nets are minimal. 

An interesting question concerns the differences be- 
tween the deterministic and random scale-free nets, even 
when the random construction is unbiased. It is evident, 
from our last example, that differences do exist. For in- 
stance, in the deterministic tree there are exactly two 
hubs (the highest degree nodes), always connected to one 
another. This need not be the case in the corresponding 
random tree. More generally, let M n (k,l) be the number 
of links with endpoints of degree 2 fc , 2 l for generation n 
of the deterministic tree. Then 



M n (k,l) =M n _i(A-l,i- 
M n (0,l)=2 l - 1 N n - 1 (l-l) 



1) 



leading to 



M n (k,l) 




0<k<l , 



k < I < n, 
k < I = n. 



(14) 



(15) 



This can be compared to M t (k,l) of random trees of an 
equivalent size (Fig. 0) . The differences seem to be triv- 
ial, leaving us with the question whether significant dis- 
crepancies show with respect to any other measure of 
structure. 
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FIG. 3: Degree-degree correlations in deterministic and ran- 
dom trees. Shown is M(4, m)/ £\ M(4, fc) as a function of 
m for deterministic trees of generation n = 8 (•), compared 
to random trees of equal size (o). The random results were 
obtained as an average over 1000 realizations. Error bars are 
smaller than the size of the symbols. 
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APPENDIX A: ON PICKING A RANDOM LINK 

How does one pick a link at random? We focus on the 
following two methods: (a) Select a link at random from 
among all the links in the network, with equal proba- 



bilities, or (b) Select a node at random, then select a 
neighbor of that node at random, and pick the link con- 
necting the two nodes. The two methods select specific 
links with generically different probabilities. 

Consider a particular link, in a graph of order N and 
size M, whose endpoints are nodes of degree fci and fc 2 . 
In method (a) each of the M links in the graph is selected 
with equal probability, so the link under consideration is 
selected with probability 



Plink 



1 

M 



(Al) 



In method (b), the probability to pick the node of de- 
gree ki first, is 1/N. Then the probability that the right 
neighbor is selected, is 1/fci. Likewise, the probability to 
select the node of degree k^ first is 1/N, and the remain- 
ing node is selected with probability l/fe. It follows that 
the probability to hit upon the specified link by random 
selection of the nodes at its endpoints is 



1 ( 1 



N \ki 



1 



(A2) 



Since 2M/N — (k), the average degree of nodes in the 
graph, it follows that 



Plink 



= W 

Pnodes 2 



1 



1 

k~2 



(A3) 



Note that this relation is valid for all connected graphs, 
regardless of structural details. An example relevant to 
our paper is a tree graph, where M = N — 1, and (k) ~ 2 
if the tree is large. In this case we deduce from (|A3() 
that only links with k\ = k^ = 2 are equally likely to be 
picked by either method. Method (b) favors links with 
min{fei, fc2j > 2. On the other hand, all the tree leaves 
(links leading to a node of degree one) are selected more 
frequently by method (a). 
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